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Abstract. We derive the number density evolution of massive field galaxies in the 
redshift range 0.4 < z < 1.2 using the A"-band selected field galaxy sample from the 
Munich Near-IR Cluster Survey (MUNICS). We rely on spectroscopically calibrated 
photometric redshifts to determine distances and absolute magnitudes in the rest- 
frame K-band. To assign mass-to- light ratios, we use two different approaches. First, 
we use an approach which maximizes the stellar mass for any _K"-band luminosity at any 
redshift. We take the mass-to-light ratio, M/ Lk, of a Simple Stellar Population (SSP) 
which is as old as the universe at the galaxy's redshift as a likely upper limit. Second, 
we assign each galaxy a mass-to-light ratio by fitting the galaxy's colours against a 
grid of composite stellar population models and taking their M/ ' Lk- We compute the 
number density of galaxies more massive than 2 x 10 10 /i _2 .M©, 5 x 10 10 /i~ 2 .M©, and 
1 x W 11 h~ 2 Mq, finding that the integrated stellar mass function is roughly constant 
for the lowest mass limit and that it decreases with redshift by a factor of ~ 3 and by a 
factor of ~ 6 for the two higher mass limits, respectively. This finding is in qualitative 
agreement with models of hierarchical galaxy formation, which predict that the number 
density of ~ M* objects is fairly constant while it decreases faster for more massive 
systems over the redshift range our data probe. 



1 Introduction 

The traditional observables used to characterise galaxies are unsuitable for study- 
ing the assembly history of galaxies, one of the most fundamental predictions 
of CDM models, since these observables may be transient. The best observable 
for this aim is, in principle, total mass, which is on the other hand very hard to 
measure. It has been argued that the best available surrogate accessible to di- 
rect observation is the near-IR ZGband luminosity of a galaxy which reflects the 
mass of the underlying stellar population and is least sensitive to bursts of star 



formation and dust extinction (Rlx fc Rieke, 1993; Kauffmann fc Chariot, 1998 



Brinchmann & Ellis, 2000). The main uncertainty involved in the conversion of 
iGband light to mass is due to the age of the population, amounting to only a 
factor of two in mass uncertainty for populations older than ~ 3 Gyr. 

The galaxy sample used here is a subsample of the MUNICS survey, selected 
for best photometric homogeneity, good seeing, and similar depth. Furthermore, 
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in each of the remaining survey patches, areas close to the image borders in 
any passband, areas around bright stars, and regions suffering from blooming 
are excluded. The subsample covers 0.27 square degrees in V (23.5), R (23.5), I 
(22.5), J (21.5), and K (19.5); the magnitudes are in the Vega system and refer 
to 50% completeness for point sources. 

The final catalog covers an area of 997.7 square arc minutes and contains 5132 
galaxies. The fields included in th is analysis are S2, S3f5-8, S5, S6, and S7f5-8. 
See Table 1 in Drory et al. (2001b ) for nomenclature and further information on 
the survey fields. 

The distances to the galaxies are derived using spectroscopically calibrated 
photometric rcdshifts. A comparison of spectroscopic and photometric redshifts 
is shown in Drory et al. (2001a ). 

We derive stellar masses by converting rest-frame if-band luminosities to 
mass using two different approaches to model the mass-to-light ratios of the 
galaxies. We discuss the resulting integrated stellar mass functions at different 
mass limits and their evolution with redshift. 

We assume J?m = 0.3, Qa — 0.7 throughout this work. We write Hubble's 
Constant as Hq — 100 h km s _1 Mpc -1 , using h = 0.65 unless the quantities in 
question can be written in a form explicitly depending on h. 



2 The maximum PLE model 

The integrated stellar mass function n(M. > -Mii m ), the comoving number den- 
sity of objects having stellar mass exceeding A4y lm , is computed using the Vmax 



formalism as desctibed in Drory et al. (2001a) 



To compute the stellar mass of a galaxy, we first use an approach which 
maximises the stellar mass for any i^-band luminosity at any redshift. 

Noting that M. / Lk is a monotonically rising function of age for Simple Stellar 
Populations (SSPs), we find that the likely upper limit for M /Lk is the mass-to- 
light ratio of a SSP which is as old as the universe at the galaxy's redshift. This 
is the most extreme case of passive luminosity evolution (PLE) one can adopt. It 
corresponds to a situation where all massive galaxies would be of either elliptical, 
SO, or Sa type. 



We take the mass-to-light ratios from the SSP models published by Maraston 



(1998 ), us ing a Salpeter IM F. S imilar dependencies on ag e are obtained from the 
models of Worthey (1994 ) and Bruzual fc Chariot (1993| ) although the absolute 



values of M./Lk vary somewhat, partly due to differences in the models them- 
selves but mostly due to the way stellar remnants are treated by the different 
authors. 



The resulting integrated mass functions for M\i m = 2 x 10 10 /i 2 -Mq, M. 



lira 



5 x lO w h~ 2 M Q , and Mum = 1 X lO n h~ 2 M0 are shown in Fig. || along with 
the integrated luminosity functions for comparison. 

The mean values of M/ Lk in the maximum PLE model in the four redshift 
bins are 0.99, 0.83, 0.73, and 0.65, as computed from the look-back time in 
our cosmology. With these mean values the mass limits correspond to absolute 
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if-band magnitudes of -22.43, -22.63, -22.77, and -22.90, respectively, for 
Mnm = 2 x 10 10 h- 2 M Q . For M\ im = 5 x 10 w h- 2 M & the numbers are -23.42, 
-23.62, -23.76, and -23.89. Finally, for Mn m = 1 x lQ lx hT 2 M® we have 
-24.18, -24.38, -24.51, and -24.64 (magnitudes with respect to ft, = 1). 

The upper and middle panels of Fig. ^| compare the evolution of the inte- 
grated luminosity to the integrated mass. It is evident that the number density 
of luminous if-band selected galaxies does not evolve significantly (given our 
uncertainties) to z = 1.2. However, because of the inevitable evolution of the 
mass-to-light ratio with z, the number density of massive systems does change. 
Transforming luminosities into masses with our maximum PLE scheme yields 
a roughly constant number density for our lowest mass limit, 2 x 10 10 h~ 2 A4q, 
and a decrease of the number density with redshift by a factor of ~ 3 for a 
mass limit of 5 x 10 10 h~ 2 \Mq, and by a factor ~ 6 for objects more massive 
than 1 x 10 n /i~ 2 .M©. As the true M/Lk at high redshift will most likely be 
lower than in our maximum PLE model, the true number densities are likely to 
decrease more rapidly with redshift. 

The steepening of the curves with increasing limiting mass in the maximum 
PLE models (despite them all having the same mass-to-light ratios at any given 
redshift) is due to the invariance of the LF with redshift and its steepness at 
the bright end. At increasing limiting mass, one is moving down the steepening 
bright end of the LF, so that the same change in the mass-to-light ratio yields a 
higher change in the number density. 

To investigate the effect the uncertainties in the photometric redshifts have 
on the values of the integrated mass function, we have performed Monte-Carlo 
simulations. The errors of the mean values of the integrated mass function (size of 
open symbols in Fig. ||) are derived by repeating the mass function analysis using 
subsamples of the template SED library, as deficiencies in the templates are the 
main source of concern for the accuracy of the redshifts during the photometric 
redshift determination. 



3 The fitted mass-to-light-ratio model 

To obtain a more realistic estimate of M/Lk, we used our VRIJK color infor- 
mation and the photometric redshift to fit the age and SFR of each galaxy using 
a grid of composite stellar populations (CSP) with 9 exponential star formation 
timescales, r, ranging from 0.1 to 10 Gyr with spectra extracted for 28 ages, 
t, between 0.04 Gyr and 15 Gyr for each r. The input SSPs for constructing 



the composite stellar populations models are taken from Maraston (2002 ), again 
using Salpeter IMF. 

Fig. [| shows the evolution of the if-band mass-to-light ratios as a function 
of age for each value of r in the grid. Except for the two largest values of r, the 
slope at ages t > 2 Gyr is remarkably independent of the actual star formation 
timescale. Moreover, the slope of the time evolution of M./Lk is the same even 
for the shortest value of r, 0.1 Gyr, which essentially represents an SSP. Let 
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aside normalisation effects, we therefore may expect a similar result for the 
mass function as obtained with the PLE model. 



3.0 




age [Gyr] z 

Fig. 1. Evolution of the K-band mass-to-light ratio, M/Lk, with age for differ- 
ent star formation histories (left panel; see text). The average isT-band mass-to- light 
r&tio,M/ Lk , of the MUNICS sample as a function of redshift as determined from fit- 
ting CSP models to the V, R, I, J, K colour data base. The dashed line denotes the 
mass-to-light ration of the maximum PLE model (see text; right panel). 



The average X-band mass-to-light ratio of the galaxy population determined 
by applying this fitting procedure is shown in Fig. [I]. The figure also shows the 
PLE mass-to- light ratio as a function of z. Apart from the different normalisa- 
tion, the evolution with redshift is very similar, a consequence of the insensitivity 
of Ai /Lk on the star formation history. 

Finally, the lower panel of Fig.|| shows the integrated mass function for the 
same mass limits as those applied above, using the individually fitted M/Lk 
values. 

The most striking feature of Fig. || is the similarity of the maximum PLE 
and the CSP-fitted curves. Note that there is a difference in the normalisation of 
the two, and due to the log scaling of the figure, the slope appears to be different 
in the plot. 

If star formation played an important role at z ~ 1, the presence of young 
populations would have pushed M./Lk down, and therefore the CSP-fitted 
curves would be expected to be steeper than the maximum PLE curve, which 
assumes no star formation happens at all after z = oo. 

Nevertheless, the number density of massive systems seems to decline, with 
this decline being stronger for more massive systems, and therefore one is inclined 
to think that merging does play an important role. Indeed, Le Fevre et al. (2000| ) 
derive a number of 0.6 to 1.8 major mergers per L* galaxy since z ~ 1 from HST- 
based pair counts of galaxies with known redshifts selected from the CFRS. We 
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Fig. 2. Comoving number density of objects having rest-frame K-band luminosities 
exceeding -22.42 + 51og/i (2 x 10 10 /i" 2 L Q ), -23.42 + 5 log ft (2 x 10 10 /i" 2 L Q ), and 
-24.17+5 log ft (2xlO lo ft~ 2 L ) (upper panels) and comoving number density of objects 
having stellar masses exceeding A4i im = 2 x 10 10 ft _2 X Q , Mi\ m = 5 x lO lo ft _2 7W , and 
-Miim = 1 x lO 11 ft _2 A / l0 (integrated stellar mass functions; middle panels). Mass to 
light ratios are assigned to maximise the stellar mass at a given luminosity (see text), 
and thus are likely upper limits. The lower panel shows the integrated mass function 
for the same mass limits, this time individually determining M/Lk for each object 
by fitting against a grid of CSP models (see text). The solid points denote the values 
measured separately in each survey field, the open circles denote the mean values over 
the whole survey area. The size of the open circles is chosen to represent our estimate 
of the total uncertainty in the mean values. 



observe a decline in the number density by a factor of > 2 for somewhat more 
massive systems, and almost no significant density evolution at L*. 

Therefore, we are inclined to think that if merging is the dominant factor in 
increasing the mass of these if-selected massive galaxies, most of the merging 
has to be dissipationless, involving rather low star formation activity. 

The main uncertainty in these conclusions is still the field to field variation, 
in spite of the relatively large area surveyed, followed by the choice of SED 
templates used in the photometric redshift code (see above). The size of the 
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open symbols in Fig. |2| represents our estimate of the total uncertainty of the 
mean values. If we assume a Gould IMF ( Gould ct al., 1998| ) instead of a Salpeter 



IMF, the evolving M/Lk curve becomes lower in its normalization as the mass- 
to-light ratio becomes smaller due to the reduced number of low-mass stars. The 
slope does not change significantly. 

The observed density evolution as a function of mass is qualitatively con- 
sistent with the expectation from hierarchical galaxy formation models. Most 
rapid evolution is predicted for the number density of the most massive g alaxies 



while the nu mber density of L*-galaxies should evolve much less. E.g. Baugh 



et al. (1998| ) predict that the number density of galaxies of a stellar mass of 
10 10 h~ 1 A4o decreases by a factor of ~ 3.1 over redshift range 0.4 < z < 1.2 (for 
the cosmological parameters as used here). Though this agreement is encourag- 
ing, both more elaborated models and improved sets of data are required. The 
latter can be obtained by larger and deeper samples, and more realistic estimates 
of M./Lk based on spectroscopic observations of the galaxies.. 
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